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Abstract: We study the one-loop level renormalization of the electrodynamics of spin 1/2 
fermions in the Poincare projector formalism, in arbitrary covariant gauge and including 
fermion self-interactions, which are dimension four operators in this framework. We show 
that the model is renormalizable for arbitrary values of the tree level gyromagnetic factor 
g within the validity region of the perturbative expansion, ag 2 < 1. In the absence of 
tree level fermion self- interactions, we recover the pure QED of second order fermions, 
which is renormalizable only for g = ±2. Turning off the electromagnetic interaction we 
obtain a renormalizable Nambu-Jona-Lasinio-like model with second order fermions in four 
space-time dimensions. 
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1 Introduction 

Second order spin 1/2 fermions were considered by Feynman for the first time in [1], fol- 
lowing a seminal work by V. Fock [2]. The V-A structure of weak interactions proposed by 
Feynman and Gell-Mann was motivated by the existence of chiral irreducible representa- 
tions of the Lorentz group for spin 1/2 fermions that naturally obey a second order equation 
of motion [3]. Feynman-Gell-Mann formalism is specially useful in the world-line formu- 
lation of perturbative quantum field theory (see [4] for a review and further references). 
The relativistic quantum mechanics and the quantum field theory of the Feynman-Gell- 
Mann equation were studied in [5-16]. The non-Abelian version was considered in [17]. In 
[18, 19], second order fermions were implemented in the lattice. Recent developments on 
the formalism can be found in [20, 21]. At one-loop level there are some partial results in 
[11, 13, 16, 18]. 

In [22], an alternative second order formalism for spin 1/2 fermions was presented, 
based on the projection onto irreducible representations of the Poincare group [23-26], an 
idea that follows from previous attempts to solve the ancient problems of the quantum 



point function ffj 
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description of interacting high spin fields. In this formalism, there exists a deep connection 
of the gyromagnetic factor of spin 3/2 fields with their causal propagation in an electro- 
magnetic background [26], and with the unitarity of the Compton scattering amplitude in 
the forward direction [27]. When applied to spin 1 fields, a similar connection between 
unitarity of Compton scattering in the forward direction and the gyromagnetic factor is 
found. Besides, in that case the gyromagnetic factor is found to be intrinsically related to 
the the electric quadrupole moment of the field [28]. Furthermore, it was recently shown 
that in this formalism the multipole moments of particles with spin 1/2, 1 and 3/2 trans- 
forming in different representations of the Lorentz group, are fixed by the tree level value 
of the gyromagnetic factor [29] and the electric charge. 

As shown in [22], the second order formalism for spin 1/2 fermions yields the same 
results as Dirac QED for the tree level amplitude of Compton scattering, whenever g = 2 
and the parity conserving solutions for the external fermions are used. However, unlike 
the spin 1 and 3/2 cases, in the case of spin 1/2 the value of g is not constrained by 
unitarity arguments. This makes the formalism interesting in the formulation of effective 
field theories for the electromagnetic interactions of hadrons where the low energy constants 
are precisely the free parameters in the Lagrangian. Thus, the study of the renormalization 
properties of this formulation of QED with arbitrary values of the gyromagnetic factor is 
of paramount importance. 

Beyond the tree level approximation, a first approach to the one-loop structure of the 
second order QED of 1/2 fermions in the Poincare projector formalism was presented in 
[30]. From the analysis of the superficial degree of divergence, only the 2-, 3- and 4-point 
vertex functions are superficially divergent. The complete 2- and 3-point vertex functions 
and the divergent piece of the only 4-point vertex function appearing at tree level (.f/77) 
where calculated in the Feynman gauge and shown to be renormalizable at one-loop level for 
arbitrary values of g. In this work we generalize the calculations done in [30] to an arbitrary 
gauge and complete the analysis of the renormalizability of the model, calculating in an 
arbitrary gauge the divergent piece of the remaining 4-point vertex functions generated at 
one-loop level. Furthermore, since the second order fermion fields have mass dimension 1 in 
four space-time dimensions, in addition to the interactions generated by minimal coupling 
arising from the gauge principle, point-like four-fermion interactions are also dimension-four 
and gauge-invariant operators. In this work we also take these potentially renormalizable 
Nambu-Jona-Lasinio-like terms [31, 32] into account and study the one-loop structure of the 
second order electrodynamics of spin 1/2 self-interacting fermions in the Poincare projector 
formalism and in an arbitrary covariant gauge. 

This paper is organized as follows: In Section 2, we present the Feynman rules and 
Ward-Takahashi identities used in the paper. In Section 3 we calculate the renormalized 
2- and 3-point vertex functions, and the divergent piece of all the 4-point vertex functions, 
closing with the derivation of the beta functions of the theory. We discuss our results in 
Section 4 and our conclusions are given in Section 5. Conventions, useful identities and 
some technical details of the calculations are presented in appendix A. Finally, the scalar 
functions arising in the calculation of the three point vertex function, at one-loop level, are 
given in appendix B. 
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2 Feynman rules and Ward-Takahashi identities 



Fermion fields have dimension 1 ((cZ — 2)/2 in d space-time dimensions) in the second order 
formalism based on the Poincare projectors [22, 30], thus local U(l) gauge symmetry and 
naive renormalization criteria allows for Nambu-Jona-Lasinio-like fermion self-interactions 
in the Lagrangian, which in this formalism are dimension-four operators. The most general 
dimension four U(l) gauge invariant Lagrangian in this framework is 



\f^F^ + D^T^D^ - m 2 ^ + ^ (^) 2 + y (VW) (M 



(2.1) 



where = 8^ + ieA^ (fermion charge — e) is the covariant derivative, Xj (j = 1,2,3) are 
the couplings of the three possible dimension four fermion self-interaction terms, and the 
space-time tensor T^ v is given by 

T »u _ gl xv _ igM nv _ ig 'M^, (2.2) 

with M» v = e^ ua/3 M a ^/2 (see [30] and appendix A for further conventions). In eq. (2.2), 
g' parameterizes parity violating electric dipole interactions, and g can be identified with 
the gyromagnetic factor. In the following, we restrict our analysis to parity conserving 
interactions setting g' = 0. 

Fixing the gauge through the non-gauge invariant contribution 

-^-gauge — 2£ ' \^-°) 

allows us to write the Feynman rules shown in figure 1 (with explicit Latin spinor indices 
running form 1 to 4). For closed fermion loops, the familiar factor of —1 given by Fermi 
statistics must be included. Here l a j, is the Kronecker delta function for spinor indices. 

Concerning fermion self-interactions, in order to have control of the different fermionic 
label combinations, we adopt an unconventional graphical device, sometimes used in the 
literature for this purpose (see e.g. [33]), which requires some comments. We unify the 
three kinds of self- interactions into a single diagram, depicted in figure 1. We remark 
that the dashed line in that diagram does not correspond to a particle exchange, it is 
only a notational convention that indicates the order in which the tensor product of the 
currents is done, and therefore, the diagram should not be seen as a reducible one [33]. 
The advantage of this notation is the automatic remotion of ambiguities when dealing with 
the contraction of spinor indices in diagrams involving self- interactions. For example, the 
self-interaction contribution to the tree level four fermion vertex function is given simply 
by the standard anti-symmetrization prescription for the fermion lines in the diagrams of 
figure 2, and yields i{\ a bcd ~ Kbad)- 

In this work we study the renormalizability of the model as it is, and we do not 
include an additional 1/2 factor in closed fermion loops. This 1/2 factor is used to provide 
a connection between the second order formalism and Dirac theory [20, 21], and we will 
see below that it is only justified in the case g = ±2, Xj = 0. 
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b P a 



iS«b(p) = i-^ 




Figure 1. Feynman rules for the second order QED of self- interacting fcrmions for an arbitrary 
covariant gauge £. 

p h b p[,a pi,b p[,a 



Pi- 



Pi- 



P2, 



Pi- 



Figure 2. Feynman diagrams for the self-interaction contribution to the tree level four fermion 
vertex function. 

Gauge invariance impose relations among different Green functions. The relevant 
Ward-Takahashi identities were derived in [30]. The first one is 

q^^p,q,p + q) = S'^ip + q) - S'^ip). (2.4) 

where —ieT^^p + q,p,q) is the jf/7 vertex function and iS'~ 1 (p) stands for the inverse of 
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the fermion propagator. We will use also this relation in its differential form 

r^o,p) = ^M. (2.5) 

The second Ward-Takahashi identity relates the jf/77 vertex to the ff"f vertex 

q^^u(p,q,p',q') = r v (p + q,q',p') - ^uip,q',p' - <?), (2.6) 

or in differential form 

r (p v' q 1 ) - dT ^ p > q '> p,) + dT ^v'>p') (2 7 ) 

L,,v(p,u,p,q)- ^ + Q ^ . [Z.<) 
The tree level vertices in figure 1 satisfy these relations. 

3 Renormalization 

In [30], the superficial degree of divergence for this theory was studied -in the absence 
of self-interactions- concluding that only vertex functions with at most four external legs 
can be ultraviolet divergent. That conclusion does not change when we introduce fermion 
self- interactions. Thus, the proof that the theory is renormalizable requires to work out all 
vertex functions up to four external legs. We will carry out this analysis at one-loop level 
for arbitrary covariant gauge in dimensional regularization, using the naive prescription 
for the chiral operator 7 s (commuting with M^ v in d dimensions, see appendix A). In 
this work, 7 s is present only in diagrams with self-interactions, and therefore, pure QED 
diagrams are free from possible dimensional regularization inconsistencies. 

3.1 Counterterms 

The parameters of the bare Lagrangian are the fermion mass mo, the fermion charge eo, the 
gyromagnetic factor go and the self-interaction couplings Aoj (j = 1, 2, 3). The renormalized 
fields are related to the bare ones as 

A? = Z~^A%, Vv = zf^ Q . (3.1) 

It is convenient to split the Lagrangian into its free and interacting parts 

jS? = JS?o + -S$ (3.2) 



where 



with 



J% = --Fg v F^ v - -L(^A) M ) 2 + d^d^ - mgVWo, (3.3) 
4 2£ 

^int = -ieo^oTo^^Vo - d*V 7Wo]^o + e^oV^Av 



J*» = gV» _ ig Q M^. (3.4) 
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Writing the free Lagrangian in terms of the renormalized fields we get 



Jgfo = - ~ A F? v F rilv - —{d»A ril f - -F^F rilv b x (3.5) 



1 „.„._ 1 ,„., . ,o 1 



Here we used the following definitions: 



5± = Z\ — 1, (^2 = ^2 — 1, #m = Z m — Z2, Z m = — ^-^2, (3.6) 



lilt 



and £ r = Z x x ^o- Similarly, the interacting Lagrangian can be rewritten as 

= -ieSrT^d^r - d^ r T r ^ r ]A v r + e 2 r ^ r A^A rfl (3.7) 

+ ~Y {tPrlpr) 2 + (VV7 5 Vv) 2 + (iprM^lpr) 2 

-ieSrT^d^r - d^ r T r ^l/) r }A»6 e 

-ie r [ipri-igrM^d^r - d^ri-igrM^r] A u r S g + e 2 r ^ r A»A rp ,h 

H 5 (VvVv) H 5 (VV7 Vv) H ^ {A M ^A) ) 



where 



and 



- 1, ^3 = ^3 - 1, <Jaj = ^ - 1 <5 9 = Z e9 - Z e , (3.8) 



Z e = —ZfZz, Z3 = -|Zl^2, ^Aj = Y^-^2> -^eg = — Z e . (3.9) 

The renormalized space-time tensor If is defined in terms of the renormalized constant 
g r as 

T» v = g^ - ig r M^ '. (3.10) 
In d = 4 — 2e dimensions, the renormalized parameters must be modified as follows: e r — > 

Notice that we have adopted a slightly different notation to the one used in [30] for 
the above definitions, in order to make our results easier to compare to Dirac and Scalar 
QED. Here the perturbative expansion is given around e r = 0, e r g r = and X r j = 0. The 
allowed values of g r by the perturbative expansion are discussed in section 4. 

In the following, for the sake of clarity, we will drop the suffix r in the renormalized 
parameters, keeping the suffix for the bare quantities. In this notation, the Feynman 
rules for the renormalized fields are given in figure 1, and the corresponding rules for the 
counterterms are shown in figure 3. 

3.2 Vacuum polarization 

The vacuum polarization at one-loop level was calculated in [30] in Feynman gauge and 
in the absence of the fermion self-interacting terms. These calculations are not modified 
by the introduction of self-interactions nor by the consideration of an arbitrary covariant 
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b P a 



i(p 2 - m^&ita,, - im 2 5 m l ab 



H 1 v 



q,H 




p,b 
q.p 



p.b 

Pub 



p'.a 

4-y 




P ,a 
P\,a 



ieK(P,p')8 e + igK b V (p'-p)v5 s ] 



p 2 ,d p' 2 ,c 



Figure 3. Feynman rules for the counterterms in the second order QED of self-interacting fermions. 




3,V q,n 




q. V 



Figure 4. Feynman diagrams for the vacuum polarization in the QED of second order self- 
interacting fermions at one-loop. 



gauge. For the sake of completeness here we just quote the relevant results. From figures 
1 and 3, the vacuum polarization is given by 



iIF"(g) = -Hl*»(q) - id 1 [g^q 2 - q»q v ) , 



(3.11) 



where — iH*^ u (q) stands for the contribution of the one-loop diagrams in Figure 4. These 
diagrams yield the polarization tensor 



U*^(q) = (g^q 2 - q^q^ir* (q 2 ) 



(3.12) 
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with 



e 2 r (3g 



24vr 2 



2m 



B {q 2 , m 2 ,m 2 ) + — [B (q 2 , m 2 , m 2 ) - B o (0, m 2 ,m 2 )] 



1 



3 

(3.13) 

Here r = Tr[l] = 4. We will systematically leave this trace unevaluated until the very end. 
It will prove to be very useful to use r as a parameter to compare our findings with well 
established results below in Section 4. 

We remark that in this work we do not attempt to correct this function with an 
additional 1/2 factor as in [30], where this was the only vertex function in which the extra 
1/2 factor was used. 

The Passarino-Veltman scalar integral Bq in eq. (3.13) is defined as 



(2ttm) 



A-d 



ITT' 



<fl 



1 



[I 2 — to 2 ] [(Z + p) 2 — to 



(3.14) 



Using d = 4 — 2e and the conventional Feynman parameterization, this function can be 
written as 



B (p 2 ,mim 2 ) = ^ + B (p',mi,mi), 



1 



2 „.2 2> 



with 



and 



1 1 



B (p 2 ,mlm 2 ) 



dx In 



o 



7 + In 4-7T 



m 2 (l — x) + m\x — p 2 x(\ — x) 
J 2 



(3.15) 
(3.16) 

(3.17) 



For future purposes, our conventions for the Passarino-Veltman scalar integrals Co and Do 
are 



ri / 2 2 2 2 2^ 



(27T/Z) 



4-d 



dH 



[I 2 - m 2 ][(l + Pl ) 2 - m 2 ][(l + p 2 ) 2 - m 2 



21 ' 



r-,/222 2 2 2 2\ 
(27T/i) 4 - d 



cfZ- 



[Z 2 - to 2 ] [(Z + pi) 2 - to 2 ] [(I + P2 ) 2 - to 2 ] [(Z + p 3 ) 2 - m 2 ] ' 
From eq.(3.11), the total vacuum polarization tensor is given by 



n^(q) = (g^ q 2 -q^Mq 



with 



<q 2 ) =Tr*(q 2 ) + S 1 . 
Using the on-shell renormalization scheme requires the form factor to satisfy 

vr(g 2 -+ 0) = 0, 



(3.18) 
(3.19) 

(3.20) 
(3.21) 

(3.22) 
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which in turn fixes the value of the counterterm as 



*i = -7r*(g 2 = 0) 



(4tt) 2 

Then, the physical form factor is given by 



I 9 1 - - 
4 3 



1 m 
- — In— = 

e (i 2 



21 



Aq 2 



e 2 r ( /3g 2 - 4 2m' 



+ — [Bo(q 2 ,m 2 ,m 2 ) - B Q (0, m 2 , m 2 )] 



24vr 2 

3.3 Fermion self-energy 

The fermion self-energy at one-loop level reads 

- i^abiP 2 ) = -i^lbiP 2 ) + i-iP 2 ~ m 2 )5 2 t ab - im 2 5 m t ab , 



(3.23) 



(3.24) 



(3.25) 



where — iY>* ab (p 2 ) is computed with the one-loop diagrams contained in figure 5. The on- 



(3.26) 



J ab\ 

shell renormalization conditions for this Green function require the propagator 

S lP) = ~ 2 \ i ■ 

to have a simple pole at p 2 = m 2 , which impose the following renormalization conditions 

0£(p) 



T,(p 2 = m 2 ) = 0, 



dp 2 



0. 



These relations fix the counterterms in eq.(3.25) as 

SV=m 2 ) 

Om — n , 02 — 



dH*(p 2 ) 



dp 2 



and the renormalized fermion self-energy is given by 



- *£(p 2 ) = -i [E*(p 2 ) - £* (m 2 )] + t (p 2 - m 2 ) 
The calculation of diagrams 1-4 in figure 5 yields 



(3.27) 



(3.28) 



(3.29) 



16vr 2 



2 {p 2 + m 2 ) B (p 2 ,m 2 ,m 2 ) + *m 2 J3 (0, m 2 , m 2 ) + ^j^m 2 



+ (1 " (P 2 - ™ 2 ) C (0,/; m 2 , m 2 , m 2 ) + (1 - £) (p 2 - m 2 ) B o (0, m 2 , m 2 ) 



m 

16vr 2 



[( r _ l)Ai - A 2 - 3A 3 ] J3 (0, m 2 , ?n 2 ) + (r - l)Ai - A 2 + ^ |> , (3.30) 



and therefore, the counterterms in eq.(3.28) read 



(4vr) 2 



1 to 2 
- - m — 

e /x z 



(r - l)Ai - A 2 - 3A 3 - 3 ( 1 + 9 — ) 



A, 



+ (r-l)Ai-A 2 + ^- 7 + ^ e 



9 2 \ „2 



(3.31) 
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p,b < l + p S p,a 




p,a 



p,b ,' l + p \ p,a 




p,a 



Figure 5. Feynman diagrams for the fermion self-energy at one-loop. 






p,b p',a p,b p',a p,b p',a 

1 2 3 



q.p 




Figure 6. Feynman diagrams for the ff^ vertex at one-loop. 



(4vr)' 



r(3-0 



1 , m 2 , m± 
m — = — m — tt 



(3.32) 



The renormalization constant of the fermion field Z2 does not depend on the gyromag- 
netic factor nor on the self-interaction couplings. Here and in the following we use a small 
photon mass m 7 to regulate infrared divergences. Finally, from eqs. (3.25,3.31,3.32) we get 
the renormalized fermion self-energy as 



S(P 2 



(47T) 5 



2(p 2 + m 2 ) [B (p 2 , m 2 , m 2 ) - B {m 2 , m 2 , m 2 )] + 4(p 2 - m 2 ) 



(3.33) 



777/ 

+ 2{p 2 - m 2 ) In — £ + (1 - £) (p 2 - m 2 ) z C (0,p 2 ; m 2 m 2 m 2 ) 



2^2 
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3.4 //7 vertex 



From the Feynman rules of figures 1 and 3, the ffj vertex at one-loop level is 

~ ieT ab(P' 1,P') = -ieV^(p,p') - ieT* a %(p, q,p') - ieV^ b {p,p')5 e - ie[igMj£ 'q v ]8 g , (3.34) 

where —ieF^(p,q,p') stands for the contribution from the one- loop diagrams in figure 6. 
It can be shown that the one-loop contribution satisfy the Ward-Takahashi identity 



q»r*(p, q ,pi) = -x*(p' 2 ) + i;*(p 2 ). 



(3.35) 



It was pointed out in [30] that the contribution of diagrams 1,2 of figure 5 and 1-3 of figure 
6 do indeed satisfy this equation. Diagrams 4,5 of figure 6 are proportional to M^ v q v and 
vanish upon contraction with q^, while the related contribution of diagrams 3,4 of figure 
5 is constant and does not modify the right hand side of eq.(3.35). Writing eq.(3.35) in its 
differential form 

dZ*(p 2 ) 



K(pAp) 



(3.36) 



and using eqs. (2.5,3.34) we get 



5 e = 5o 



(47T) 



(3-0 



1 77i mr, 
- - m — - In — £ 



(3.37) 



From eqs. (3. 8, 3. 9), this relation also fixes the counterterm for the //77 vertex function 



£3 = 5 2 . 



(3.38) 



The total contribution of diagrams 1-5 in figure 6 can be written as 



r* M (p, q,p') =E*l<f + F*lr^ + &*igM^q y + WigM^r v 
+ I*igM a ^r a qp + TigM aP r a q p q*, 



(3.39) 



where E*-JJ* are scalar form factors and = (p + p')^. A convenient decomposition of the 
form factors is the following: 



<D* = ^OiP^, (3.40) 

i=0 

with 0* = E*,F*,G*,H*,F, J*. Here (z = 0, 9) are scalar functions and PVi denote 
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the following Passarino-Veltman scalar integrals: 



PV 
PVt 
PV 2 
PV 3 
PVi 
PV 5 
PV 6 
PV 7 
PV 8 
PV 9 



n/n 2 /2 2 2 2 2\ 

D (0,p ,p ;m 7 ,m 7 ,m ,m ), 

r< (n 2 2 2 2\ 

C (0,p ;m 7 ,m 7 ,m ), 
C (0,p /2 ;TO 2 ,m 2 ,m 2 ), 

/-i / 2 /2 2 2 2\ 

Co(p ,p ;m 7 ,m ,m ), 

B (q 2 ,m 2 : m 2 ), 

B (p 2 ,m 2 ,m 2 ), 

B {p >2 ,m 2 ,m 2 ), 

B (0,m 2 ,m 2 ), 

B (0,m 2 ,m 2 ). 



The explicit form of the scalar functions (Dj is presented in appendix B. 

The form factors E*,H*,F and JT* turn out to be finite. Furthermore, E*,H*,JJ* 



vanish on-shell (p 2 = p 



12 



m 



q 2 = (p' — p) 2 = m 2 — > 0). The remaining form factors 



take the following values in this limit: 



F 



OS 



GT 



OS 



(4vr) s 



(3-0 



(4vr) 



1 , m 2 mi 

- - In — - In — 77 

e fi z m z 

1 m 2 

- - In — =- 



(3.41) 



l-^)e 2 + A 1 + A 2 



1 + i> 



2 

(3.42) 



I, 



OS 



2(4vr) 2 m 



2 ' 



(3.43) 



There is a misprint in [30] for I os which does not affect any other result. The on-shell 
renormalized vertex function in eq.(3.34) reads 



-ie (l + 5 e + F os ) Ir" - ie (1 + d e + 6 g + G os ) i 5 M^^ + I os ^M^r a ^r^. 

(3.44) 

The counterterm in eq. (3.37) cancels the divergence of the charge form factor and yields 
Fos = 1 f° r the corresponding on-shell renormalized form factor. This choice for 5 e also 
cancels one of the divergences of the magnetic form factor. In fact, the coefficient of the 
egM^ v q v term in eq.(3.44) reads 



1 + 5 e + 5 g + G* os 
1 



1 + Sa + 



(47T) 



1 m* 
- — In— ^ 

e p, z 



l-^)e 2 - 



1 + 



+ 2e 2 + ^ 



(3.45) 



We choose the following value for the 5„ counterterm to remove the remaining divergence 



1 



(47T) 



1 

- — In— ~ 
e /z z 



l- 9 -)e 2 ^ 



1 + 



(3.46) 
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According to this choice, the on-shell value of the renormalized magnetic moment form 
factor is given by 

a_ A 3 
2vr 32vr 2 



9Gos=g[l + ^ + ^), (3.47) 



with a = e 2 /(4-7r). 

In summary, the renormalized vertex function in eq.(3.34) reads 

V =Etq> M + Elr^ + GigM^q v + MigM^ u r u 
+ teg M a ^r a q^ + Mg M ap r a q^, 



(3.48) 



with the finite form factors 

E = E*, H = H*, 1 = 1*, I = T, (3.49) 
and the renormalized form factors 

F = 1 + F*-F5 S , (3.50) 

G = 1 + £ + 3^ +G *- G °s- (3 ' 51) 

3.5 //77 vertex 

The //77 vertex function at one-loop level is obtained from the Feynman rules of figures 
1 and 3 as 

ze 2 0, q,P>, q') = ie 2 V% + ^e 2 ^^(p, q,p>, q>) + ie 2 V^8 3 , (3.52) 

where the one-loop corrections ie 2 r*^ (p,q,p' ,q') are given by the diagrams of figure 7. 
The counterterm 5s has been already fixed in eqs.(3.32,3.38). 

In [30], it was pointed out that the contribution of diagrams 1-3 of figure 6 and 1-9 of 
figure 7 satisfy 

q^(p,q,p',q') = T*> '(p + q,q> \p') - T* v (p,q',p' - q). (3.53) 

It can be easily shown that this relation is unmodified with the inclusion of the remaining 
diagrams in figures 6 and 7. Indeed, diagrams 4,5 of figure 6 depend only on = (p' — p)^ 
and their contribution to the right hand side of eq.(2.6) vanishes. On the other hand, 
the contribution of diagrams 10-15 in figure 7 satisfy q^T*^ v {p, q,p', q')\ 1015 = 0. There- 
fore, eq.(3.53) holds for the full set of one-loop diagrams in figures 6 and 7. Using now 
eqs. (3.38,3.53) in eqs. (3.34,3.52) it can be explicitly shown that the second Ward-Takahashi 
identity in eq.(2.6) holds for the renormalized vertex functions. 

The divergent pieces of the loop contributions to the r*^(p, q,p', q') vertex function 
can be isolated taking the zero external momentum limit. In this limit, the sum of the first 
two diagrams can be written as 

ze 2 rr(0,0,0,0)| 1+2 = -^2^[e + ^]il ab + O( e °). (3.54) 
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q',v q,li 




p,b 



10 



q'.V q.fl 




p.h 



,,.v 




p,b p' \a 
12 



q-li q',v 

p + q + l 
p,b \p + l p' + l' p',a 



<f,v 



p+i'^-^^y+i 

p.b > v p'-q + l i p'.a 
14 



q,H 



p,b \p + l p' + l; p'.a 



Figure 7. Feynman diagrams for the ffyy vertex at one-loop. 



Similarly, the divergent piece of the third diagram is 



ie J r:f(o,o,o,o)| 



2g^[t + ^-}-l ab + 0(e»). 



(3.55) 



(4tt) 2 ^ L '' 4 J e 

Notice that this divergence cancels the one coming from the first two diagrams in eq. (3.54), 
yielding a finite contribution of the first three diagrams. The calculation of the next pure 
QED contributions is straightforward: 



ie 2 rX(0,0,0,0)| 



ie 



(4vr) : 



;2{?, + i)g^-t ab + 0{et 



I 6 " 9 =™4£ s f»li ab + O(e ) 
(47r) z e 



ie'VJf (0,0,0,0)| 



(3.56) 



(3.57) 



Adding up eqs. (3.54,3.55,3.56,3.57), we obtain the divergent part of the pure QED loop 
divergent contributions to the 77// vertex function as 

ic a rr(o,o ( o,o)| 1 - 9 = -^ 5 '- 1 



(3-e);[2ieVlla b + 0(e u ). 



(3.58) 
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There is no room for further contributions to 83 due to the fermion self-interactions. 
Therefore, the renormalizability of the model requires the finiteness of diagrams involving 
self-interactions in figure 7. For these diagrams one has 



ab 



• 2-n*V I 1 l\ 1 10—12 2\ 2<e 

ie T£ (p,q,P,q)\ =e KbcdV 



d d l f[V v (l,l-q')V»(l-q,l)\ 



dc 



+ 



(2vr) d i □[/-?]□[/]□[/ -g' 
lV»(lJ + q)V»(l + q',l)] dc 2g^l dc 



ie T ab (p,q,p>q 



1 ,,113-15 



U[l + q'}U[l}U[l + q] 
d d l 



c 2 \ ,, 2<E 



a[i-q f ]D[ 



q) 



(3.59) 



+ 



(2vr) d i D[l + p]D[l + p'} 
[Vn£ - q + l,p' + l)V(p + l,p'-q + Ok 

D\p + l]D\p' + l]D\ P '-q + l] 
[V(p + q + l,p' + l)V»{p + l,p + q + l)] dc 



(3.60) 



U[p + l\U[p' + l}U[p + q + l] 
Here and in the following, we adopt the shorthand notation □ [p] = p 2 — m 2 and A [q] 



Again, the possibly divergent part of this contribution can be isolated taking 



vanishing external momenta and is given by 
10-15 _ o f d d l 



ie 2 T*^ 



ab (0,0,0,0) 
= 2eW e 



2eV e 



d d l 
(2vr) c 



(27r) d 
1 



AM 1 



^-dc (^abed ~ ^cbad) 



[t - ljAi - A 2 A 3 



(3.61) 



1 



ab ■ 



This integral is finite. As a consequence, the sum of the diagrams involving fermion self- 
interactions in figure 7 is free of divergences. Thus, the divergent part of the full set of 
one-loop diagrams in figure 7 is 

,,2 



ie 2 rX(0, 0, 0, 0) = -— ^ (3 - ? [2ieVll afe + <V). (3.62) 



This divergent piece is proportional to V^i and is exactly of the same magnitude and 
opposite sign to the divergent piece in £3 in eq.(3.38), which has been already fixed from 
the Ward-Takahashi identities. Inserting eq.(3.38) in eq.(3.52) we get a renormalized .f/77 
vertex function which is free of ultraviolet divergences. 

The calculations so far presented generalize the results obtained in [30] to an arbitrary 
covariant gauge and the inclusion of fermion self-interactions. It is easy to show that all 
the above results reduce to those of [30] in Feynman gauge £ = 1 and vanishing fermion 
self-interactions Xj = (J = 1,2,3). Now we turn our attention to the 4-point vertex 
functions not considered in [30]. 

3.6 47 vertex function 

The four gamma vertex function is absent at tree level and is generated at one-loop level 
by the diagrams in figure 8. Notice that the fermion self-interactions do not contribute to 
these diagrams and there is no counterterm for this vertex function. Thus, if the model 
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Figure 8. Fcynman diagrams for the 47 vertex at one-loop. There are 9 additional diagrams 
obtained from diagrams 1-9 reversing the charge flow in the loop. We denote these diagrams with 
a prime, e.g., V stands for diagram 1 with arrows pointing in the opposite direction 



is renormalizable, the sum of all these pure QED diagrams must be finite. Again we are 
only interested in the study of the renormalizability of the theory and we will focus on the 
possible divergent parts of these diagrams. We can isolate the divergent pieces setting all 
the external momenta to zero. In this limit, the triangle diagrams 1 — 6 and 1' — 6' of figure 
8 (see figure caption for an explanation of the primed notation used here) yield 



-ie r^ Q/3 (o, 0,0,0) 



1-6,1' -6' 



8re 4 -U 2e 
a 



d d l 
(2ir) c 



n[if 



rpfiva/3 



with 



The divergent piece of the box diagrams 7 — 9 and 7' — 9' of figure 8 is 



-ie 2 r^ a/3 (o, 0,0,0) 



7-9, 7'-9' 



-Are" 



24 



,2c 



d(d + 2 



■I'' 



and for the remaining diagrams we have 

10-12 



- ; e 



! r^ Q/3 (o, 0,0,0) 



-4reV e 



d d l 
{2ir) d 

d d l 



(/ 



2\2 



rpfiua/3 



(2vr) c 



□ [J] 



rpfivaf} 



(3.63) 
(3.64) 

(3.65) 
(3.66) 
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pi,b p[,a pi,b p[,a p\,b p[,a p\,b p[,a 




P2,d P2,c P2,d p' 2 ,c P2,d p' 2 ,c P2,d p' 2 ,c 

13 14 15 16 



Figure 9. Feynman diagrams for the 4/ vertex function «A* bcd (pi,P2,Pi,P2) & t one-loop. There 
are 16 additional diagrams obtained from the crossing {p[,a} O {p'2,c}, which contribute to 

i A *cbad(Pl,P2,P2,Pl)- 



Adding up the contributions in eqs. (3.63,3.65,3.66), all the divergences cancel. Therefore, 
the four-gamma vertex function is finite, as expected. 

3.7 4/ vertex function 

The last potentially divergent vertex (according to the analysis of superficial degree of 
divergence), is the four-fermion vertex function (three photon vertex function vanishes 
because of charge conjugation). From figures 1 and 3, this function is iA a b c d (pi,P2,Pi,P2)~ 
i^cbad (Pi , P2 , p'i 1 Pi) > witn 

iKbcd {pi,P2,p'l,P2) =iKbcd + iKbcd {PUP2,Pl,P2) ,„ c7) 

+ i [5 Xl Xttabtcd + 5 X2 \2ll b ll d + 5 Xs A 3 MjM Ml/cd ] . 

Here «A* fecd is obtained from the loop diagrams in figure 9. 

We use the same technique to isolate the possible divergent terms in these diagrams. 
We start with the pure QED diagrams. In the zero external momenta limit, diagrams 1,2 
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of figure 9 yield 



iA: 6cd (0,0,0,0)| 1+2 = 2e 



,2 , (d - l)g 4 



/' 



(2tt) c 



2\2 



Similarly, the divergent piece of the triangle diagrams 3,4 is 



3+4 = _ 4e 4 



*A^(0, 0,0,0) 
Finally, for diagram 5 we get 



,2 , (d - l)g 2 
£ H ^ 



.■lc 



(2tt) c 



□ [/] 2 A[/] 2 



/ 2 

□[J]A[J]s 



(3.68) 



(3.69) 



iA: fccd (0,0,0,l 



^)4e 4 [e 2 + (d-l)] ^/ 



(2tt) £ 



A[i]J 



(3.70) 



Notice an additional symmetry factor 1/2 in this last contribution, the same that appears 
in Scalar QED. This factor is not related to the extra 1/2 factor in fermion loops used 
in [20, 21] (and should not be confused with it). Instead, this 1/2 is the appropriate 
factor required to correct the double-counting of internal photon configurations that comes 
from the factor 2 contained in each V^ v vertex. In appendix A, we give a more detailed 
calculation of the above results. 

The divergent pieces of the integrals in eqs. (3.68,3.69,3.70) are alike. Integrating and 
adding up these contributions we get 



iA* abcd (0, 0,0,0)| 



1-5 



• 4 

ie 



(4vr) 2 



2\ 2 i 



6 1-^ ^t ab l cd + 0(e»). 



(3.71) 



In addition to the 47 vertex function, which turns out to be finite and dictated by pure 
electrodynamics (no influence of the fermion self- interactions), the 4/ vertex function must 
be finite in order to conclude that the pure QED studied in [30] is a renormalizable theory. 
eq.(3.71) is a gauge independent result, and besides, is no affected by the 7 s issues related to 
dimensional regularization. Thus, from eq.(3.71), we conclude that pure electrodynamics 
of second order fermions is a renormalizable theory only in the case g = ±2. As we 
will discuss in the next section, these values of g lie in the range where the perturbative 
expansion is valid. This sharp prediction for the electrodynamics of second order fermions 
in the Poincare projector formalism is modified by the inclusion of fermion self-interactions, 
which relax this condition and provide the adequate counterterms to have a renormalizable 
theory for arbitrary values of g. This is an analogous situation to the case of scalar QED, 
where a similar divergence appear at one loop level in the 4-scalar vertex function. The 
quartic — ^X^^ <p) 2 coupling is necessary in scalar QED to remove this divergent dynamical 
interaction. 

So far, all the renormalized vertex functions are finite at one-loop level for self- 
interacting fermions. The proof of the renormalizability of the QED of self- interacting 
fermions requires the calculation of the diagrams involving fermion self-interactions in fig- 
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ure9. The corresponding divergent pieces of these diagrams are 

,2 



iA: 6cd (0,0,0,0)| 6 - n =e 2 <|2eA abcd + ^- 



X a ecfM^M^ fd + X ebfd M^M i 



fiucf 



+ XefcdM^M^fb + XabefMgM^fd 

ie 2 1 



k 



I 2 



d d l 



(2Tr) d \A[i]n[i} 2 



(An) 2 e 



\g 2 (Ai + A 3 ) + 2^A X 



2<? 2 (A 2 + A3) + 2£A 2 



tabled 



+ 



<? 2 (2Ai + 2A 2 - A 3 ) + 2£A 3 



2 

12-16 



M^M^ + O(e ) 



iA: fecd (0,0,0,0)| 

{ — XabefXfed + X aec fX eb fd + X ae fdX ebc f + X ae f b X e f c d + A a fe e j X ce fd) ^ J 

= J^l\{ [( 4 " ^) A ? + 2A^ + 3A§ + 2AiA 2 + 6A!A 3 ] Wed 
+ [(2 - r)A 2 + 3A 2 + 6AxA 2 + 6A 2 A 3 ] 7^ 



+ 



--^A! + 6AiA 3 + 6A 2 A 3 



} + 0(i 



(3.72) 



d d l 



[2n) d \D[l] 



(3.73) 



Adding up all contributions, we get the total divergent piece of the four-fermion vertex 
function as 



iA* abcd (0, 0,0,0) 



i 1 



(4tt) 2 e 



2\ 2 



1-2- e 4 + e 



3o 2 

-|- (Ai + A 3 ) + 2£Ai 



+ (4 - r)A 2 + 2A 2 + 3A 2 + 2A X A 2 + 6A1A3 }t ab t c d 



i 1 



(4^) 2 e 
i 1 
(Z^j 2 ? 

+ O(e ). 



+ 



3^ 
2 

„2 



(A 2 + A 3 ) + 2£A 2 



+ (2 - T )A^ + 3X 2 + 6AiA 2 + 6A 2 A 3 



(2Ai + 2A 2 - A 3 ) + 2£A 3 



4 + T 



A3 + 6A1A3 + 6A 2 A 3 > Mj%M, 



L ab 1V1 ^cd 



(3.74) 



From eq.(3.67) it is clear that, in the case of self-interacting fermions, we still have at our 
disposal the 5\ . counterterms to absorb this divergence. The choice 



1 1 



(4vr) 2 e 



3g 2 



l + ^)+2£ 



.Ai . „A 



(3.75) 



+ (4 - r)Ai + 2-4 + 3-^ + 2A 2 + 6A3 
Ai Ai 



- 19 - 



1 1 



J \2 



? A 3 



(4tt) 2 e 

1 1 

~1 



\2 

+ 3-^ + 6Ai + (2 - r)A 2 + 6A 3 
A2 



(4vr) 2 e 



1+2$ 



+ 6A1 + 6A 2 



4 + r 



A; 



(3.76) 



(3.77) 



yields a finite four-fermion vertex function and completes our analysis the renormalizability 
of the QED of second order self-interacting fermions at one-loop level. 

3.8 Beta Functions 

Summarizing, form the results obtained in eqs. (3. 23, 3. 31, 3. 32, 3. 37, 3. 38, 3. 46, 3. 75, 3. 76, 3. 77) 
and the definition of the counterterms in eqs. (3.6,3.8,3.9), the relation between the bare 
and renormalized parameters of the theory is given by 



e = Z x 2 Z 2 Z e fi e e, 



Z-y 1 Z 2 1 Z 3 n Ae e z , Xoj — Z 2 "ZxM 2e X 



,2e„2 



2c ' 



'3 ' 



5o — Z e 1 Z eg g, 



ITln 



(3.78) 



with the renormalization constants (defined all the MS subtraction scheme, for consistency) 



MS 



Z 



<7 MS 
£ 2 



7 MS 



e 2 r / g 2 



(4vr) 



r,MS _ ^MS 



-)- 

3) r 

e 2 1 
1 + — (3-£)r, 



(4vr) 5 



+ 

2 V Ai 



+ 2$ 



(3.79) 
(3.80) 
(3.81) 



A 2 A 2 II 

+(4 - r)Ai + 2-^ + 3-^ + 2A 2 + 6A3 }-, 
Ai Ai e 



7 MS 
J A 2 



(4vr) 



(3.82) 



A 1 

+3-^ + 6A1 + (2 - r)A 2 + 6A3 }-, 
A2 e 



7 MS 
J A 3 



7 MS 



rMS 



1 



1 



(4vr) 2 1 
+6A1 + 6A 2 

rMS _|_ ^MS 



r f Ai a 2 

^ ( as + A3 



1+2$ 



1 + 



1 



(4tt) 



2-e+f 1 



-Ai - A 2 + ( 1 + - ) A 3 



r^MS 1 cMS 

^5 + m 



1 + 



(47T) 



(r - l)Ai - A 2 - 3A 3 - ( t\ + ^- ) e 2 



3^ 



(3.83) 



(3.84) 



-. (3.85) 
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According to these constants, the two different relations between eo and e in eq.(3.78) 
collapse to 

e = Z~ l/2 fi e e, (3.86) 

just as in Dirac and Scalar QED. 

From eqs. (3.78-3.85) one can extract the following beta functions (3 n = and 
anomalous dimensions 7 m = in the e — > limit: 

e 3 r /3 



48tt 2 \4 



.9 



32^ 2 



.2 /„2 



9 — 4) — 4(Ai + A 2 ) + 4 ( 1 + - ) A 



_L/2 

16vr 2 | 4 



e 4 ( 92 _ 4) 2 +3e2[(4 + 92) Ai+5 2 Aa 



(3.87) 
(3.88) 

(3.89) 



/3a 2 



£a 3 



In 



+2(4 - r)A 2 + 4A 2 (Ai + A 2 ) + 6A 3 (2Ai + A 3 ) 

- I ^{3e 2 [(4 + 5 2 ) A 2 + 5 2 A 3 ] 

+12A 2 A! + 2(2 - t)X\ + 6A 3 (2A 2 + A 3 ) }, 
--±^{e 2 [(l2-g 2 ) A 3 + 2 5 2 (A 1 + A 2 )] 
+12A 3 (A 1 + A 2 )-(4 + r)A|}, 

{-3e 2 {g 2 + 4) + 4[(r - l)Ai - A 2 - 3A 3 ]} . 



64vr 2 



(3.90) 



(3.91) 



(3.92) 



As expected, all £ dependence drops out and the beta functions are gauge invariant. Finally, 
taking r = Tr[l] = 4 in eqs. (3.87-3.92), we obtain the definitive one-loop beta functions 
and anomalous dimensions of the theory: 



e 3 /3 2 



12vr 2 V 4 
9 



32vr 2 

= ~ 16^*1 



[e 2 {g 2 - 4) - 4(Ai + A 2 - 3A 3 )] , 
1 ^-/(g 2 -4) 2 + 3e 2 [(4 + g 2 )X 1+g 2 X, 



(3.93) 
(3.94) 

(3.95) 



Px 2 
/3a 3 

7m : 



+4A 2 (Ai + A 2 ) + 6A 3 (2Ai + A 3 ) 
1 



64vr 2 



2 {3e 2 [(4 + g 2 ) A 2 + g 2 X 3 ] + 4A 2 (3Ai - A 2 ) + 6A 3 (2A 2 + A 3 )} , (3.96) 
^ {e 2 [(12 - g 2 ) A 3 + 2g 2 (A x + A 2 )] + 4A 3 [3(Ai + A 2 ) - 2A 3 ]} , (3.97) 
{ -3e 2 {g 2 + 4) + 4[3Ai - A 2 - 3A 3 ] } . (3.98) 



1 



Notice that f3 g in eq. (3.94) vanishes for {g = ±2, Ai = A 2 = A 3 = 0}, and for g = with 
arbitrary Xj . The analysis of the evolution of the coupling constants is beyond the scope 
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of this paper and here we will just relate these special fixed points of j3 g to the well known 
Dirac and Scalar QED, and the simplest model of pure self interactions. 



4 Discussion 

As stated in the previous section, (3 g vanishes for g = ±2, Ai = A2 = A3 = 0. In this case, 
the theory has a simple connection to Dirac QED, as expected. It amounts to make the 
replacements r — > (l/2)Tr[l] = 2, g — > ±2 and Xj — > in eq. (3.87) to obtain the well 
known Dirac-QED beta function and the corresponding anomalous dimension of the mass 



Here the additional 1/2 factor in r is the same used in [20, 21], which can be traced to the 
relation between the Dirac operator and the Poincare projector operator for g = 2: 



The recovery of Dirac QED shows that the perturbative expansion for g is justified in the 
g = ±2 case, and therefore, for —2 < g < 2. This observation is important because g 
is a running coupling constant in this theory. One could be tempted to assume that the 
perturbative expansion is done around e = and g = 0, however, as stated in Section 
3.1, the true couplings in the Lagrangian are e and eg, thus, the perturbative expansion is 
given by powers of e 2 /(47r) = a, and (eg) 2 /(47r) = ag 2 . This double expansion may lead 
to terms of order ag, as -for example- the one-loop corrections to the gyromagnetic factor 
Ag = ag /2tt in eq. (3.47). In general, the validity of the perturbative expansion will be 
driven by the conditions a<l, ag 2 <C 1. Hence, even if we consider values like e.g. g = 5, 
at low energies we will have a perturbative expansion in powers of ag 2 ~ 25/137 = 0.18, 
which is still a reliable parameter. Notice, however, that for g ^ ±2 the Xj terms are 
generated at one loop level. Said in other words, if we put Xj = at some energy scale 
(say atomic scale, for example), the pure QED dynamics generates self-interactions at 
one- loop level at a different energy scale, as shown in eqs. (3.93-3.98). Therefore, the 
evolution of g itself is influenced by its implicit self-interaction coupling dependence, and 
this dependence must be taken into account even in explicitly Xj independent results, like 
the case of f3 e . In this theory, all couplings e, g and Xj run with the energy in an intricate 
pattern, and the naive estimation of the perturbative regime must be modified. These 
modifications are expected to be small for small values of Xj at low energies, but certainly, 
the elucidation of the energy scale where the perturbative expansion is valid requires to 
solve the renormalization group equations. 

Similarly, we should be able to recover results of scalar QED under the appropriate 
considerations. The vanishing of (3 g at g = points to this possibility, and it is compatible 
with the fact that for this special value of g, there is no spin dynamics and the spin degrees 
of freedom must collapse to a multiplicative factor. This is indeed the case whenever we 
also change the Fermi statistics of our fermion degrees of freedom to the Bose statistics of 
conventional scalars. Formally, the theory described by eq. (2.1) can be turned into scalar 



D 
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QED if one changes the statistics and rescales the degrees of freedom and couplings accord- 
ingly. At the beta-function-level, this is accomplished with the following replacements in 
eq. (3.87): g — > 0, A2, A3 — > 0, r — > — 1 (Fermi— >-Bose and reduction of degrees of freedom) 
and Ai — > — ^ (rescaling and sign rectification of quartic coupling in scalar QED). This 
procedure yields the desired scalar QED beta and gamma functions: 



Finally, we remark that this result does not mean that, in the case g = 0, the theory 
describes just several copies of scalar QED. Statistics play an important role here. Taking 
g = in eq. (3.93) changes completely the nature of the theory, since the sign for (3 e flips: 

^=0 = -^. (4.4) 

Thus, in this case (and in general for g 2 < 4/3), our results point to an asymptotically 
free theory. Conclusive results require to solve the renormalization group equations in eq. 
(3.93), in order to understand better the contents of the theory and the new effects produced 
by arbitrary values of g and Xj. However, it is interesting that these fixed points of j3 g -in 
the second order formalism for spin 1/2 based on the Poincare projector- correspond with 
the two simplest realizations of renormalizable theories. 

While this paper was under preparation, a work with interesting results for second order 
QED [35], was announced. In that work, a non-perturbative derivation of (3 e is given, based 
on a generalization of the minimally-coupled Klein-Gordon equation to include a Pauli term 
with arbitrary g. The j3 e derived in [35] coincides with eq. (3.93) up to a constant factor 
(related to r and the 1/2 additional factor) for —2 < g < 2. Also there, the possibility of 
asymptotic freedom behavior for fermions was recognized. 

Regarding fermion self- interactions, many results involving the A,- couplings are subject 
to the issues of 7 s in dimensional regularization. The naive prescription used in this work is 
known to be internally inconsistent, but it has proven to give the correct results in anomaly- 
free theories (see for example [34] and references therein), like ours. Previous works on the 
renormalization of two dimensional fermion models (which have deep similarities to our 
model) show that the possible problems related to the evanescent operators, that arise when 
the Clifford algebra is continued to d dimensions, start only at two-loop level [33], thus, it 
is reasonable to expect that the results involving fermion self-interactions obtained in this 
paper are reliable. Admittedly, verifications of the regularization procedure are necessary to 
draw definitive conclusions for the renormalizability of the QED of self-interacting fermions. 

Finally, another attractive feature of the theory is the case e = 0, i.e., when we switch 
off the electromagnetic interactions. In this limit, we are left with a theory of self- interacting 
fermions, similar to the one proposed long ago by Nambu and Jona-Lasinio, but which in 
our formalism turns out to be renormalizable. There is a special class of models of this kind 



- 23 - 



in which dimensional regularization also gives unambiguous results. Taking, for example, 
e = 0, Ai = A, A2 = A3 = 0, we obtain the following two-parameter model: 

jgf = d^d^ijj - m 2 ^ + - (^V) 2 , (4.5) 
with one-loop beta and gamma functions 



Px = 0, 7 m = — 2> ( 4 -6) 



3A 

16^ 

which is free from 7 s inconsistencies. Besides, the massless limit of this model is indeed 
finite in dimensional regularization at one- loop for arbitrary values of A. 

5 Summary and conclusions 

In this work, we studied the one-loop level renormalization of the electrodynamics of sec- 
ond order fermions in the Poincare projector formalism, in an arbitrary covariant gauge 
and including fermion self-interactions. In contrast to Dirac fermions, the second order 
ones have dimension mass dimension 1 in d space-time dimensions). Thus, four- 

fermion interactions are dimension-four operators and -according to the superficial degree 
of divergence- they must be included at tree level. There are three main conclusions of this 
work. First, if we start with vanishing tree level self-interactions (pure QED), at one-loop 
level all the Green functions with up to four legs turn out to be renormalizable, except 
for the four-fermion interaction generated by the loops, which is renormalizable only in 
the case g = ±2. Hence, pure QED of second order fermions in the Poincare projector 
formalism is renormalizable only for these values of g. This sharp prediction is obtained 
using dimensional regularization but it is not affected by the known inconsistencies of this 
regularization method related to the definition of 7 s in d space-time dimensions. Second, 
the introduction of fermion self-interactions at tree level modify this behavior and render a 
renormalizable electrodynamics of second order self-interacting fermions for arbitrary val- 
ues of the gyromagnetic factor. In this case, g becomes a true coupling constant, running 
with the energy. This theory may be used for composite particles like baryons, where the 
gyromagnetic factor is one of the low energy constants and its physical value is a suitable 
starting point for the formulation of the corresponding effective field theories. Third, if we 
turn off the electromagnetic interaction we are left with a four-dimensional renormalizable 
model of second order self-interacting fermions, which may be relevant in the formulation 
of effective theories in the strong regime, along the proposal by Nambu and Jona-Lasinio, 
but with second order fermions. In general, the last two main conclusions involve the 
7 5 issues of dimensional regularization method, and further verifications by alternative 
regularization methods would be desirable. In the case of the third conclusion, however, 
there is a class of models with scalar-scalar self-interactions which is also free of possible 
inconsistencies of dimensional regularization. 
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A Notation and useful identities 



Lorentz generators satisfy the following algebra in 4 dimensions: 

[ M a£ = _i^ M ^ _ g ^M^ + gP v M a » - g^M au ), (A.l) 

{M^, M af} } = hg^g^ - g^g va )t + ^ Q/3 7 5 , (A.2) 
( e 0123 = — eoi23 = 1)- The chirality operator 7 s is defined as 

7 5 = '-M^M^, = \e^M a ^ (A.3) 

and satisfies 



[ 7 5 ,M^]=0, (7 5 ) 2 = 1. (A.4) 

Higher products of the generators can be calculated using recursively these relations. 
We also need to calculate the trace of the product of generators. The simplest one is 

Tr [MH = 0, (A.5) 

as required from Lorentz covariance. Similarly, for 7 s we have 

Tr [ 7 5 ] = 0, Tr [ 7 5 M^] = 0. (A.6) 

In d dimensions we assume that the generators still satisfy the Lorentz algebra in eq. 
(A.l) and the anti-commutator relation in eq.(A.2), but now with g^ = d. Also in d 
dimensions eq. (A.5) and Tr[l] = 4 hold. The well known issues of 7 s in dimensional 
regularization come from the impossibility of satisfying simultaneously both eq.(A.4) and 



Tr 



7 5 M Q/3 M^ ^0. (A.7) 



As Tr [7 5 M a/3 M^] is nowhere needed in our calculations, we use the naive dimensional 
regularization prescription, which amounts to keep eqs.(A.4,A.6) unmodified in d dimen- 
sions. 

The following reduction formulas for tensor integrals are useful in the determination 
of divergences: 

d d l m v i a i? f d d l (I 2 ) 2 . 

(A.8) 



(27r) d A[l] m D[l] n d(d + 2)J (2vr) d A[l] m D[l] n ' 
d d l IH V f d d l I 2 



J 



(A.9) 



(2vr) d A[l] m \3[l] n d J {2-K) d A[l] m U[l} n ' 

with given by eq. (3.64). 

As eq. (3.71) constitutes a crucial result, here we work out in detail the divergent part 
of its constituent diagrams. Thus, defining P a ^(q,C) = g a ^q 2 — (1 — £) 9 \ , for diagrams 
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1,2 of figure 9, the divergent part is given by 
iA* abcd (0,0,0,0)\ 1+2 

„4..4e / d d l P w (l,£)PvAl,Z)[Vni, 0)^(0, l)] ab [VP(-l,0)V°(0,-l) + V°(l,0)Ve(0,l)] cd 



e fi 



(2vr) c 



n[l] 2 A[l]< 



d d l 2e{l 2 ) 2 t ab l cd + gXl^F (M^M^) ab {M w , M ua } cd 



eV 4e 



(2vr) a! 

d d l 2{l 2 ) 2 
(2vr) d n[/] 2 A[/] 2 



□ [/] 2 A[Z] 2 



ft ab t cd + 



2d(d + 2) 



ab 



' cd 



(A.10) 



Here eq. (A. 8) was used in the last step. Notice that this result, and in general all pure 
QED divergent parts, are unaffected by the 7 5 issues of dimensional regularization, as 
the {Mp P ,M va } contributions always appear contracted either by l p l a or by g pa , giving 
unambiguous results when continued to d dimensions. This fact is more evident when 
the generators are written in terms of the conventional Dirac matrices. The evaluation of 
the product of generators in eq.(A.lO) can be performed with the aid of eqs.(A.l,A.2) as 
follows: 



{M pa M u a ) ab {M/, M up } cd + (m^M^) ^ { {AV, M vfi } + {M^M UQ \ } cd 



(M^M\) ab 
d(d-l) 2 3d(d-l) 



l )g pv \ cd 



ab 

+ ( M^M^ 

V J ab 



ua ) 



•-cd 



+ 



taht 



ab^cd 



d(d + 2)(d- 1) 



t a bt cd . 



(A.ll) 



In this way, eq.(3.68) is obtained from eqs. (A. 11, A. 10). Similarly, for diagrams 3,4 of 
figure 9, the use of eqs.(A.l,A.2,A.9) gives 

*A: bcd (0,0,0,0)| 3+4 

d d i p w {i, OPuAi, i)g pa {\y»(i, o)V(o, i)] ab led + [v»(i, 0)^(0, i)] cd t ab } 



-2eV 4e 



-2eV 4e 



-4eV 4e 



-4eV 4e 



(2vr) d □[/]A[/] 2 
d d l 2ei 2 t ab tcd + 9 2 U P [{M^Mp P ) ab t cd + {M^M^) cd t ab ] 



d d l J 2 

(27r) d a[l]A[l] 2 

d d l I 2 
(2vr) d D[l]A[l] 



n[i]A[if 



flatted + h [(M^ a M m ) ab t cd + (M^M pa ) cd t abi 



2 1 ? ^ 2d 



9_ 

2d 

2 r 2(f(d-l) 



tabtcd, 



in agreement with eq.(3.69). Finally, the contribution of diagram 5 of figure 9 is 



iA^O, 0,0,1 



2) eV \/ 



d d l P w (.l,0PuAl,0 [2g^tab] [2g pa t cd ] 



A[l}< 



(A.12) 



(A.13) 
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and reduces to eq.(3.70) upon contraction of Lorentz indices. 

We close this appendix with a list of some useful products needed in the evaluation of 
eqs.(3.72,3.73): 

[M^M a %[M F %] d = \ ab t cd + ^abl'cd - 1M%M„ vcd + 0(e) (A.14) 

[M^M a % b [M a pM^] cd = h ab l cd + ^ilYll d + 2M%M^ cd + 0(e) 

M^M^M^ = -M aP + 0(e) 
[l b M^] ab [^M^] cd = M^M^ cd + 0{e). 

B Scalar functions for the decomposition of the form factors of the three 
point function ffq 

The scalar functions O, = Ei,Fi,Gi,Hi,Ji,Ii from the decomposition of the form factors 
in eq. (3.40) are the following: 

E =0, 

Ei =4C(1 - (P 2 ~ P' 2 ) {P 2 ~ m 2 ) (p' 2 - m 2 ) (p ■ p' + m 2 ) , 

E 2 =4C(1 - (p 2 - m 2 ) (p' 2 - m 2 ) (p ■ p' + p 2 ) , 

E 3 = - 4C(1 - {p 2 - m 2 ) [p' 2 - m 2 ) (p • p' + p' 2 ) , 

£ 4 (p 2 - p ' 2 ) | ( g 2 - 4) [m 2 q 2 +p-p' (p 2 + p' 2 ) - 2p 2 p' 2 } 

+ 8 [m 4 + 2p-p (m 2 + (p ■ p )) - p 2 p' 2 ] } 

- 4C(1 - £) {P 2 ~ P' 2 ) {P 2 ~ m 2 ) (p' 2 - m 2 ) , 
E 5 =-C{p 2 - p' 2 ) [{g 2 -A)q 2 + S(p-p' + m 2 )] , 

E 6 =C [{g 2 ~ 4) {p 2 - p' 2 ) {p 2 -p- p') + 8p 2 {p' 2 + m 2 ) + 8p ■ p' (p 2 + m 2 )] , 
E 7 =C [(g 2 - 4) (p 2 - p' 2 ) (p' 2 -p-p') - 8p' 2 (p 2 + m 2 ) - 8p ■ p' (p' 2 + m 2 )] , 

^8 =0, 
Eg =0. 
^0=0, 

Fi =4C(1 - OQ 2 {p 2 ~ m2 ) {p' 2 ~ m2 ) {P-P' + m2 ) > 

F 2 = - 4C(1 - {p 2 - m 2 ) [p ■ p' (q 2 -p 2 - m 2 ) + p 2 (p' 2 + m 2 )] , 

F 3 =- 4C(1 - £) (p' 2 - m 2 ) [p ■ p' (q 2 - p' 2 - m 2 ) + p' 2 (p 2 + m 2 )] , 

Fa =(q 2 { {g 2 ~ 4) [m 2 q 2 + p ■ p' (p 2 + p' 2 ) - 2p 2 p' 2 ] 

+ 8 [m 4 + 2p ■ p' (m 2 + (p ■ p') ) - p 2 p' 2 ] } , 

F B =- <q 2 [{g 2 -4)q 2 + 8{p-p' + m 2 )] - - £)q 2 (p 2 - m 2 ) (p' 2 - m 2 ) , 
F G =C [{g 2 - 4) {p 2 - p • p') q 2 + 8p- p' (p 2 - m 2 ) - 8p 2 (p' 2 - m 2 )] , 
F 7 =C [(g 2 - 4) (p' 2 - p ■ p') q 2 + 8p ■ p' (p' 2 - m 2 ) - 8p' 2 (p 2 - m 2 )] , 
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^8=0, 

F 9 =0. 
G 



(4vr) 2 

Gi = - 8C(1 - (p 2 - m 2 ) (p' 2 



m 



I A 2 2/2 

(p • p J — p p 



G 2 =8C(1 - (p 2 - m 2 ) (p • p') z - p 2 p' 2 



G 3 =8C(1 - (p' 2 - m 2 ) (p • p'Y - p 2 p' 2 



G 4 =2( <j 2m V + 2m 2 (p' 2 - p 2 



) 2 + 4p-p'{(m 2 +p-p') «? 2 -2 [(p-p') 2 -pV 2 ]} 
+2p' 4 (p-p' - p 2 ) + 2p 4 (p-p' - p' 2 ) + (g - 2) (m 2 + p • p') (p' 2 - p 2 ) 2 } , 



2C {2 (m 2 + p • p') g 2 + 5 (p 2 - p' 2 ) 2 + ( 5 2 + 4) 



p 2 p' 2 — (p • p 



f]} 



+ 



1 



(4tt) 2 

p 2 



-t{ 2 p~ 1 



Ai + A 2 

2 + p • p') (p 2 — p- p') + 2m 2 



2 /2 
p p 



(p-p')' 



+ 5P 2 {P 2 ~ P 2 ) {P 2 +P- P 1 ) } , 



r - 2C 



2p' (m + p ■ p') (p' — p • p') + 2m p p' — (p • p') 



G 8 =C 



+ 9P' 2 (p' 2 ~P 2 ) {p' 2 +p-p') }, 
4m 2 



/2 1 2\ / A z 2/2 
p 2 p ,2 \P +P)[ip-P) ~PP 



G 9 = -8C(l-0 [{p-pf-p 2 p 2 

H =0, 

Hi =0, 

#2=0, 
#3=0, 

H4=2(g (p 2 -p' 2 ) (?n 2 +p-p') q 2 
H 5 = -2(g (p 2 -p' 2 ) t? 2 , 
2C 



H s ="| {-2 (p 2 - m 2 ) [(p • p') z - p 2 p' 2 j + gpV (p 2 + p • p') I , 



#7 



#8 



2C 



p' 



/2 



-2 (p' 2 — m 2 ) (p • p') — p 2 p' 2 



+ W (p /2 + p-p')}, 



4£m 2 (p 2 - p' 2 ) 



/ A z 2/2 
[P-P) ~P P 



#9=0. 

/o =2Cg 2 , 
h =0, 
/ 2 =0, 

is=0, 
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c 



(p ' p'Y ~ p 2 p i2 



|3mV + 6q 2 P ■ P' [m 2 (p 2 + p' 2 ) - p 2 p' 2 ] 



+ 2 (p-p') 2 (p 2 + p' 2 - 2m 2 ) + p 2 p' 2 (p 2 + p' 2 - 8m- 



(q 2 



(p -p'Y — p 2 p' 2 

c 



[3 (p 2 + p' 2 ) (to 2 + p ■ p') - 2 (p • p) (p-p' + 3to 2 ) - 4p /2 p 21 



{3p 6 (to 2 + p-p' -p' T 



P' 



{p-p') 2 — p 2 p' 2 

A 



P 



9m 2 p ■ p' + p' 2 (5to 2 + 7p ■ p) — 6 (p ■ p') 2 — p' 4 
+p 2 |4to 2 (p ■ p') 2 + p' 2 —5m 2 p ■ p' — 2 (p ■ p') 2 + 2 (p • p') 3 j + 2m 2 (p ■ p') 3 j , 



h 



c 



p 



J 2 



{p-p') 2 — p 2 p' 2 



{3p' 6 (m 2 +p-p' -p 2 ) 



+ P 



1 1 



-9m 2 p ■ p' + p 2 (5to 2 + 7p ■ p') — 6 (p ■ p') 2 — p 4 



+p' 2 |4m 2 (p ■ p') 2 + p 2 — 5m 2 p-p' — 2(p-p') 2 + 2 (p • p') 3 j + 2to 2 (p • p') 3 1 , 



2 ^J[(P 2 +P 2 )P-P-2P 2 P' 2 ] 



h =o. 

J =2C ( P 2 -p' 2 ) 
Ji =0, 

h =o, 
«/ 3 =o, 



C(p 2 -p' 2 " 



J A 



■h 



{p ■ p'Y — p 2 p' 2 



{ig {p-p' 



i 2\ \( A 2 2 /2 
+ TO J (p • P J — P P 



+ g 2 (6to 2 p • p' + 3to 4 + p /2 p 2 ) 



+8to" 



(p • p') 2 - p' V + 2 (p • p') 2 (p 2 + p /2 ) - 4p /2 p 2 p • p'} 



C(p 2 -p /2 ) 

{p-p') 2 — p 2 p' 2 



{-2g[(p.p') 2 -p 2 p' 2 
3p • p' (p 2 +p' 2 )+ 4p' 2 p 2 + 2 (p • p') 2 } , 



3m 2 q 2 



c 



p 



{p ■ p'Y — p 2 p' 2 



{ 2gp 2 (p 2 + p • p') (p • p') 2 - p V 2 



-p 



3p 4 — p 2 p' 2 — 2 (p ■ p') 2 (p /2 — to 2 ) —p-p 5p' 2 p 2 — 3p 4 — 2 (p ■ p') 2 (p 2 — ?n 2 ) 



c 



p 



/2 



o '2 / /2 , A / A 2 2/2 

-25P (P +P-PJ (P-Pj -PP 



+P 



{p • p'Y — p 2 p /2 

o 14 12 2 n ( A 2] / 2 2\ , / ["i- 2 12 o /4 o / A 2 1 / 

dp — p p — lyp ■ p ) [p — m j + p • p 5pp — 3p — lyp ■ p ) [p 



12 



21 / /2 2 
TO 



2(TO / 2 I2\ 1 A 
^2"(P "P ) (P-P), 



Jfl =0. 
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Here, the global factor £ stands for 

c = — e2 



1287T 2 



(p ■ p'Y — p 2 p' 2 
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